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Abstract 



(N 

>, 

C^ , We study new supergravity solutions related to large- A'^c A/" = 1 supersymmetric gauge field the- 

^^ ■ ories with a large number Nf of massive flavors. We use a recently proposed framework based on 

^>—^ ' configurations with N^ color D5 branes and a distribution of Nf flavor D5 branes, governed by a 

function NfS{r). Although the system admits many solutions, under plausible physical assumptions 
f~^ ' the relevant solution is uniquely determined for each value of a; = Nf/Nc. In the IR region, the solu- 

c7^ I tion smoothly approaches the deformed Maldacena-Niinez solution. In the UV region it approaches 

a linear dilaton solution. For x < 2 the gauge coupling Pg function computed holographically is 
negative deflnite, in the UV approaching the NSVZ /3 function with anomalous dimension 70 = —1/2 
(approaching -3/(327r^)(2Afc - 7V/)g^)), and with /3g -^ -00 in the IR. For x ^ 2, Pg has a UV 
fixed point at strong coupling, suggesting the existence of an IR fixed point at a lower value of the 
^\ ' coupling. We argue that the solutions with x > 2 describe a "Seiberg dual" picture where Nf — 2Nc 

*~^ , fiips sign. 
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1 Introduction 

An important area of application of AdS/CFT dualities concerns J\f = I supersymmetric SU{Nc) Yang- 
Mills theory with an arbitrary number Nf of fundamental flavors. The Nf = case was constructed 
by Maldacena and Niiiiez (MN) 1 , building up on a geometry previously found in [2]. Massless flavors 
in the fundamental representation of the SU{Nc) gauge theory can be incorporated following the idea 
of |3] by adding Nf spacetime filling branes. The resulting holographic models |31 [51 [B] have led to 
many interesting physical insights, including, for instance, aspects of Seiberg duality (see also [71IH1I1]). 
However, the presence of a singularity in the IR region limits the applicability of this geometry. Recently, 
a new Af = I supersymmetric geometry has been found by Conde, Gaillard and Ramallo 1101 which 
includes the previous ones as particular cases, but more generally can circumvent the IR singularity. The 
solution is parametrized in terms of a function S{r) that encodes how fiavor branes are distributed in 
the space. The plan of this paper is to use this framework to construct new solutions and investigate 
new physical properties, with the aim of understanding the extent to which these geometries can describe 
aspects of TV = 1 supersymmetric SU{Nc) Yang-Mills theory with Nf fundamental massive fiavors. 

Af — 1 supersymmetric SU{Nc) Yang-Mills theories with Nf massless flavors, analogously to their 
non-supersymmetric counterpart, are likely to abandon the QCD-like confined phase for sufficiently large 
number of flavors and develop a conformal phase before the loss of asymptotic freedom. The restoration 
of conformal symmetry and the presence of a so called conformal window in the number of flavors would 
thus identify a new family of non abelian gauge theories which is worth to explore. From the more 
phenomenological point of view, conformal symmetry might play a relevant role in the description of 
particle dynamics at energies above the electroweak symmetry breaking scale. 

Until now, the emergence of conformal symmetry in theories without supersymmetry has been discussed in 
the context of Schwinger-Dyson equations for chiral symmetry in the ladder approximation |llj , truncated 
non-perturbative RG flows [T^], supersymmetry inspired conjectures |T3] and deformation theory [T3]. 
The proof of existence of a conformal window however depends on our ability to describe these theories 
in a non-perturbative manner, following the evolution of parameters all the way from strong coupling to 
weak coupling. Lattice studies are currently the only ones to provide a fully non perturbative analysis, 
and in TVc = 3 QCD they have recently produced evidence that Nf = 12 is plausibly close to the end- 
point of a conformal window [THl HSl IZl HH] ■ Similar results have been found in [TH] using the worldline 



formalism. 

A further insight comes from supersymmetric gauge field theories. The renormalization group of A/" = 1 
supersymmetric QCD (SQCD) has been extensively studied, and the perturbative j3 function for the 
gauge coupling is given by the well known NSVZ formula [5^ 

where 70 is the anomalous dimension of the Nf (massless) fundamental superfields Q^ and Qu- The 
formula (II. ip is unchanged, with a common anomalous dimension for all fundamental fields, as long as 
the superpotential preserves a Zjq, flavor symmetry. It also applies in the massive flavor case in the far 
ultraviolet, whereas it is drastically modified below the mass scale of the Nf flavors. The anomalous 
dimension 70 is a non-trivial function of Nf, Nc, the gauge coupling and all other couplings appearing in 
the superpotential. All the information about the dynamics of the strongly or weakly interacting theory 
is thus contained in 70. Its form depends on the precise choice of the superpotential and determines the 
presence of fixed points in the parameter space of the theory. Any rigorous prediction for the existence 
and width of a conformal window would thus require to derive the /3 functions and anomalous dimensions 
of the theory in a non-perturbative way. Holographic techniques may allow to study the renormalization 
group flow beyond perturbation theory. In particular, for the supersymmetric pure Yang-Mills theory 
(SYM), i.e. SQCD with Nf = 0, the gravity dual was constructed in [T and reproduced some interesting 
features oi J\f = 1 SYM, including the correct structure of the NSVZ /3 function (jl.l[) at Nf = for the 
gauge coupling |21[|22) . The holographic /3 function contains, in addition, non-perturbative contributions, 
and the gauge coupling goes to infinity in the IR - a realization of ordinary confinement. 

The holographic approach has some well known limitations that we will not attempt to resolve in this 
paper. One limitation is that holographic models typically include an infinite tower of Kaluza-Klein 
states which have no counterpart in the field theory that one would like to describe, e.g. super QCD 
or super YM, and should therefore be decoupled. It was argued in [51 [S] that integrating out Kaluza- 
Klein excitations leads to effective couplings - e.g. quartic couplings - in the superpotential, so that 
the final effective field theory will be an Af — 1 supersymmetric QCD with the addition of an effective 
superpotential containing higher dimensional couplings. Such couplings can lead to the disappearance of 
the conformal window. An example is the role of four-fermion operators in non-supersymmetric gauge 
theories. The presence of such operators has to be traced back to chiral dynamics and the breaking of 
chiral symmetry, with consequent disappearance of conformalitylJ 

Despite these limitations of the holographic approach, it remains important to investigate the holographic 
predictions for the gauge coupling /3-function, since they provide detailed hints on the possible structure 
of the effective field theory in the strong coupling regime. Such hints may be helpful when studying the 
theory on the lattice or by means of other techniques. In particular, for Nf — 2Nc, a prediction arising 
from our study is the existence of a non-trivial UV fixed point at some strong coupling g^. Consistency 
with the RG evolution at weak coupling requires the existence of an IR fixed point at gl < g*, as we shall 
discuss. Notice that the presence of a UV fixed point at strong coupling in addition to an IR fixed point 



^ In particular, the Schwinger-Dyson gap equation in the ladder approximation says that the onset of chiral symmetry 
breaking occurs for a critical gauge coupling where the anomalous dimension of the fermion mass operator ■jm = 1 and thus 
the four-fermion operator becomes relevant in the RG sense. In 1231 has been also suggested that, while the value of the 
critical coupling will be affected by higher order effects in the perturbative (ladder) expansion, the scaling of the fermion 
propagator with ■jm = 1 is the non-perturbative signal of chiral symmetry breaking. 



has been conjectured already in the pioneering work by Banks-Zaks [53], and might lead to a mechanism 
of disappearance of the conformal window via the annihilation of a pair of fixed points as suggested in 

This paper is organized as follows. In section 2 we review the type IIB string backgrounds recently 
constructed in [lOj . The solutions are parametrized by a profile function S{r) which for particular choices 
reproduces previous solutions in the literature. In section 3, we will adopt a convenient choice of S(r) 
and obtain numerical solutions for the three cases: i) Nf < 2Nc, ii) Nf — 2Nc and iii) Nf > 2Nc- The 
gauge coupling /3-functions of the corresponding gauge theory for all cases are derived in section 4, with 
focus on the emergence of fixed points in the RG flow. We conclude in section 5, comparing our results 
with physical expectations and discussing possible future directions. 

2 The string/supergravity background 

We consider a type IIB supergravity background with Nf D5 branes extended on a non-compact two-cycle 
of a CY3-fold. The ansatz is given by the following (Einstein frame) metric and RR three-form [TO] 

dxl^ + e^^^^'-^dr^ -f e2M'-)(rf02 + ^in^ Od^'') 

g2fc(r) 

((cji + a{r)def + {Co2 - a{r) sin 0dip)'^) H j—i^s +cosedipy 

- {uii + b{r)d9) A {Cj2 - b{r) s'mOdip) A (el's -I- cos 9d(p) 

+ (6'(r) + Li{r))dr A {-dO ACoi + sin 6 dtp A C02) + (1 - b{rf + L2{r))s\a9d9 A dip hCj^ 

Nr 

- -ysm9d9Ad(pACj3, (2.1) 

where the uji are the left-invariant forms of SU(2) given by 

uJi = cos^jdO + simp siii9d(p , 

UJ2 = — sin 'ijjd9 + cos ijj sin 9d(p , 

UJ3 = dip + cos9dip , (2.2) 

and a, 6, /, k,g, h are functions of the radius r. This background embeds a flavor deformation by means 
of a stack of flavor D5 branes parametrized by the Li{r),i — 1,2 functions. It is convenient to introduce 
the functions P, Q and the profile function S 

Q{r) = e2»(acosh(2r) - 1) , P{r) = ae^^ sinh(2r) , S{r) = - — L2(r) . (2.3) 

N j 

The inverse relations for the metric functions g{r) and a{r) read 

e^s = P coth(2r) - Q , a = ^ , ^ ^ -— — . (2.4) 

^ ' ^' Pcosh(2r)-Qsinh(2r) ^ ' 



ds" 
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The profile S{r) characterizes the embedding of flavor D5 branes and it may be interpreted as an energy- 
scale dependence of the effective flavor number NfS{r) in the description of fundamental massive flavors 
of the dual gauge theory. The BPS equations can be solved in terms of the functions in (|2.3p as follows 

(2.7) 

2h 1 P'^-Q' 2k P' + NfS{r) 

4Pcoth(2r)-g ' 2 ' ^ ' 

(f) ph+g+k 

where $(r) is the dilaton and $o a constant. The function Q is an integral over the profile S{r) 

2Nc - NfSip) 



Q = coth(2r) 



-dp + qo 



(2.10) 



/o coth^(2p) 
with go a constant of integration. Finally, the function P{r) satisfies the "master" differential equation 



P" + N,S' + ,P' + «,S, (£^^|^ + P'.Q'.mS _,„,„„,,■) ^ „ 



(2J1) 



Once the profile S{r) is determined, a solution is obtained by first computing Q{r) in (I2.10p . then solving 
p. lip for P{r). It should be noted that regularity of the geometry (see (12.8^ ) requires 

P > IQI , P' > -NfS . (2.12) 

In the particular cases S = and S — I one finds solutions which already appeared in the literature. We 
briefly review these cases before moving to the massive case described by the function S{r). 

S — 0: Unflavored solution 

Setting S" = in (|2.1ip amounts to set Nf — 0. In addition to the regular solution of the Maldacena- 
Niinez model [I] [2] , there is a one-parameter deformation found in [4 that leads to solutions with regular 
behavior at r = 0. The infrared asymptotic of this unflavored one-parameter family of solutions has been 
explicitly written for P in [5], 

„ , Ahi f 47V2\ ^ 16hi ( 4A^2 32Ar4\ 

where h^ is the one parameter that labels each solution of the family. When h^ — 2Nc, one recovers the 
Maldacena-Niihez (MN) solution 1 . It is worth noting that the resulting function Q is the same for any 
value of ft.1, 

Q = 7Vc(2rcoth(2r) - 1) . (2.14) 



5' = 1: Massless flavors 

By setting S — 1 one obtains the solutions of [1] describing massless flavors. The asymptotic for these 
solutions was discussed in full detail in [5]. For large radius, i.e. in the ultraviolet of the dual gauge 
theory, a generic solution behaves exponentially 

P = fc e^''/^ + 0(1) , (2.15) 

where k is an integration constant. There are also special solutions with the following linearly rising large 
r asymptotic: 

P=|2iV,-iV/| r + 0(l) , Nf^2N,. (2.16) 

A further analysis is required to see if the geometry can actually be extended to r — > oo or if it meets a 
singularity before. This will be discussed below. When Nf = 2Nc, there are special solutions with the 
following asymptotic behavior: 

RN 
P = Po + 0{e-n, Po^_^, (2.17) 

with 

90 = ^^^^' c = 1 + v/9^4e + e^ , 0<e<4. (2.18) 

S{r): Massive flavors 

Following [lOj one can characterize the supersymmetric D5 brane embeddings by two algebraic equations 

Fi(z,) = 0, F2{z,)^0, (2.19) 

for the four complex variables Zi (i = 1, ... ,4) parametrizing a deformed conifold, thus satisfying 

ziZ2 - zszi = 1 (2.20) 

and related to the radial coordinate through 

4 

^|zi|^ =2cosh(2r) . (2.21) 

In particular, the choice made in [10 is given by the following embedding parametrized by two complex 
constants A and B: 

z-i = Azx , Z4 = Bz2 . (2.22) 

This equation, together with ()2.20|) and ()2.2ip determines the minimum distance r^ that this embedding 
reaches 



cosh(2.,)= V^^l-l v;^l-l . (2.23) 



\l^AB\ 

It depends on the moduli of A and B, as well as their phase. By demanding that the WZ term of the 
action of the full set of D5 branes in the ten-dimensional theory coincides with the action obtained from 
the embeddings one arrives at [TU] 

, , A /cosh 4r — cosh irg 
Sir) = ^^-^— — — ^e(r - r,) . (2.24) 

V2smh(2r) 



Notice that S{r) is continuous at r = Vq, while S'{r) diverges as S'{r) ~ (r — r,)"^/^ near r,, and 
it is thus singular. To avoid this singularity Conde, Gaillard and Ramallo [10] have proposed a brane 
setup for which the tip of the branes r^ is "smeared" , so that an average should be made over brane 
distributions with different tip positions, weighted with a density function p(rq). After performing the 
change of variables y = cosh(4r) and j/g = cosh(4rq) with y > 1, and assuming that the branes are 
distributed over the whole space < r < oo, then the profile function will be given by 



S{y)^l\yqp{yq)^^^, (2.25) 

where the measure function p{yq) satisfies the normalization condition 

dygPiy,) = 1 . (2.26) 

3 Simple solutions for massive flavors 

On the gauge field theory side, one expects that the asymptotic physics for Nf massive flavors at high 
and low energies should be as follows: a) at energies lower than the flavor mass (infrared limit) it should 
converge to the unflavored case, 5' = 0; b) at high energies (ultraviolet limit) it should converge to the 
Nf massless flavors case, S = 1. This picture can be realized by the gravity dual background when the 
function S{r) interpolates between the infrared/small radius limit S'(r) — >■ for r <^ rq - with rq being a 
measure of the common quark mass - and the ultraviolet /large radius limit S{r) — > 1 for r ^ rq. 
Thus we are interested in solutions for massive flavors that approach the deformed MN solution (|2.13|) 
in the infrared, i.e. in the small radius limit r — >■ 0. In [TU], to describe flavors with a given mass 0{yq) 
with some spread, a measure function p{yq) with a finite support around yq was chosen. Here, we slightly 
depart from this approach. Given the freedom in the choice of distribution of branes, we will conveniently 
adopt a smooth distribution p{yq) of branes, chosen to meet the following requirements: 

• S{r) is assumed to be a monotonous, continuous function varying between and 1, approaching 1 
at infinity. We demand S{r) ^ r^ (or smaller) for r ^ 0, so that the curvature invariants of the 
geometry near r = are the same as in the deformed MN solution. In this way we ensure that the 
metric is regular at the origin (and that there are no massless flavors). 

• In order to have a more tractable differential equation (|2.1ip . we demand that S is such that the 
integral (j2.10p defining Q can be explicitly performed with a simple result for Q. 

• Finally, we demand that p{yq) is positive definite and satisfies the normalization condition (j2.26p . 

We found an extremely simple choice that meets all these requirements: 

S'(r) = (tanh(2r))^" , n = 2, 3, . . . (3.1) 

This corresponds to a distribution of branes with masses concentrated around the maximum of S'{r), at 



y^max = arccoth I \l _ I , (3.2) 



which increases with n (for large n, rmax ~ l/41ogn). The spread Ar decreases with n. In order to 
determine p{yq), wc note that the integral defining S is related to an Abel Transform as follows 

2dy{^V~l S{y)) = A[p{v)] = r dy, -jk£= . (3.3) 



2« r^dy{V¥^S{y)) 



The inverse Abel Transform formula is 

It is easy to veriiy that the normalization condition p.26p is satisfied for this measure function. 
For the choice p.ip . we find 

4V2T{n+l) (y,~l)"-i 

In particular 

15(y, - 1) 
P(«=2)(2/g) = ;^,, _^ ,7 ■ (3.6) 

Next, we compute Q{r) in (12.101) . The basic integral we need is 

/ dr tanh"(2r) = — ^ ^f ^Fi[l, ^{1 + m), ^{3 + m),tanh\2r)] . (3.7) 

Jo ^("^ + IJ 

For integer 771, this reduces to simple expressions. Thus we find 
Q{r) - i(2A^,-Ar/)(2rcoth(2?-)-l) 

Nf f, V^ f t&nh''^-\2r) tanh'=-i(2r) tanh'=+"(2r) \ \ 
"^l SI ^^ r^- fe + n + 1 Jj ' ^'-'^ 

Q(„=2)(0 = ^(2Ar,-7V/)(2rcoth(2r)-l) + ^tanh2(2r) + ^tanh*(2r) . (3.9) 

Notice that we have set go = 0. The reason is that the term go coth(2r) produces a singular behavior at 
r = 0, thus violating our condition that the solution reduces to the deformed MN solution at r = 0. 

In the following section we will proceed to the analysis of solutions P{r) of the master equation (|2.1ip as 
a function of the parameter x = Nf/Nc- In all cases we will use the S{r) given by p.ip with tt, = 2 and 
hence Q given by p.9p . 

In general, the resulting differential equation p. lip admits the following boundary conditions: 

^ - {z:.%x> ■ '-"■ ("«) 

\2Nc-Nf\r, xj^2 
P « { Po + e-^'', x = 2 , r>l. (3.11) 

k e^^l'^ , any x 

For each a; ^ 2, the solution to (|2.1ip . P(r), is uniquely determined if we demand the following asymptotic 
conditions: 



a) At r ^ 0, P ^ hir, i.e. the solution reduces to the deformed MN solution with the asymptotic 
behavior given by (|2.13p . 

b) At large r, the solution has the linear behavior (|2.16p . P ^ \2Nc — Nf\ r. 

For a generic integration constant hi above some critical value, the large r asymptotic behavior is P ~ 
e^'"' ^, as discussed earlier. At the critical value of hi the solution has the linear behavior P ~ {2Nc — Nf)r, 
or constant for x — 2, and at any lower hi it meets a singularity before reaching r = oo. Hence, the 
condition of linear behavior at infinity specifies the solution uniquelyo 

In order to solve the differential equation (12. lip numerically, as mentioned above we take the brane 
distribution p.ip with n — 2, and Q given in p. 91) . This describes massive flavors with a mass around 
r « 0.5 (see (|3.2p ). determined by the maximum of S'{r), shown in fig. [1] together with S{r). 

S(r), S'(r) 




Figure 1: S{r) (solid line) and S'{r) (dashed line). The maximum of S'{r) at r « 0.5 indicates the 
characteristic mass scale of the massive flavors. 

Since we have to meet boundary conditions at zero and infinity, we employ a shooting method. This 
determines the critical hi. 

Figures [2^, b,c,d illustrate the solutions in the three cases Nf < 2Nc, Nf = 2Nc and Nf > 2Nc. 

• In the first case we take x — 7/4, for which we find 

— = 1.53218706 , X = - , (3.12) 

Nf 4 ^ ' 

and the solution is reported in fig. [2^. 

• In the special case x — 2 the solution that starts with P = hir near r = and asymptotes to a 
constant at infinity has 

-^ = 1.42475837 , x = 2 . (3.13) 



^The solutions with exponential behavior at infinity have a constant dilaton and become Ricci flat, which is not the 
expected asymptotic behavior for holographic applications. Some interesting applications of these solutions as describing 



properties of 6d field theories have nevertheless been found in [6] . 





p = 


= Po - e-=('-'-i) + 0(e~4'^) , 


Po = 


32Nc 
15 


' ^"^+ 2 ' ^^ 15 ' 



The large radius behavior is given by 

(3.14) 

with 

■^9 AT , /9T 8 Ar 

(3.15) 
where ri is a numerical constant. The solution is shown in fig. ^}p. 

• Finally, fig. [2}: shows a case with x > 2, taking in particular x = 7/3, for which we find 

— =^ 1.35890843 , x = - . (3.16) 

iV/ 3 

Note that x = 7/3 is related to a; = 7/4 (used in fig. ^j^) by a; — >■ x/{x — 1), which is produced by 
the change Nc — > Nf — Nc. We have made this choice for later comparison between theories related 
by a naive Seiberg duality transformation. We will comment on this below. 

More generally, one can determine hi as a function of x, with < x < oo, as shown in fig. [3l For a; — > 
we obtain hi/Nf — ^ oo. Indeed one can verify that hi — > 2Nf/x ~ 2Nc as Nf — > 0, recovering the MN 
boundary condition at r = for P discussed above. Furthermore, we note that for large x the critical hi 
approaches a finite asymptotic value, 

— ^ 1.72102763 , a; -> oo . (3.17) 

Nf 

The reason is that for a; ^ 1, one can scale P — >■ NfP so that the master equation (I2.11[) becomes 
independent of x, as Q becomes proportional to Nf, see p.9p . This scaling solution is shown in fig. [2Ji. 

4 Gauge coupling /^-function 

To compute the /3 function of the gauge coupling in the dual field theory we need to first identify the 
gauge coupling constant in terms of geometrical quantities. For the ansatz p.ip . this has been done in 
[1]. The gauge coupling turns out to be directly related to the P function as follows 

87r2 / p25 \ 

V = 2fe2'' + ^(a-l)2j =tanh(r)P(r) . (4.1) 

The second crucial ingredient necessary to obtain any /3 function in the dual field theory is the precise 
relation between the radial coordinate r of the supergravity background and the energy scale of the gauge 
theory. Relevant discussions on this point can be found in [2 |^ [551 IU l^- In particular, in [3T1 [3^, by 
adopting the energy scale defined by the gaugino condensate - a protected, i.e. gauge invariant operator 
with no anomalous dimensions in the gauge theory - in the Maldacena-Niiiiez model, and by dimensional 
arguments, the following relation is obtained 

3 



A 



a[T) , (4.2) 
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(b) 






Figure 2: The function P{r)/Nf solution to the master equation ()2.1ip that matches between the deformed 
Maldacena-Niihez solution (|2.13p in the infrared (r — >■ 0) and the linear behavior in the ultraviolet 
(r —> cx)). The dashed line corresponds to Q{r)/Nf {\Q{r)\/Nf in fig. c). (a) x = 7/4 (b) x — 2. (c) 
X = 7/3. (d) X = oo. 



where /i is an arbitrary renormalization scale at which the gaugino condensate is defined, and A is the 
scale dynamically generated by quantum corrections in the gauge theory. Strikingly, this relation applied 
to the unflavored Maldacena-Nufiez background leads to a gauge coupling f3 function that reproduces the 
complete perturbative NSVZ formula (|l.ip for Nf = 0. The relation (|4.2p gives rise to the UV behavior 



— ^63 



r > 1 . 



(4.3) 



In extending the relation between fi and r to models with Nf ^ niassless fiavors, one needs to consider 
a number of issues. In particular, interesting solutions exist with a = 6 = in (|2.ip . so one should seek 
for other possible definitions of the energy scale than (14.21) . As emphasized in [3] fBl , for a class of flavored 
M — 1 supersymmetric models, the same UV relation (14. 3p arises from any of the following identifications 



a{r) 



b{r) 



A 



,-20(r) 



(4.4) 
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N, 
10 r 
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Figure 3: Near r = 0, the solutions are required to behave as P ~ hir to approach the deformed MN 
solution. The figure shows the critical values of the parameter hi which are required for P to have linear 
behavior at infinity. 

The relations (|4.3p and (14. 4p can be generically written in the form 

- = F{r) , F{r) -^ e^^-^ for r ^ oo . (4.5) 

Different choices of F are analogous to the ambiguity that appears on the field theory side in the choice 
of renormalization scheme. Using (|4.ip and (|4.5p . we obtain the following expression for the /3 function: 



3F 3F 

/38^2/ 2 = — — -(9r(tanhr P) = j" (sinhr coshr P' + P) (4.6) 

t _F'cosh r 

and by knowing F and the solution P, we can now compute fi^^iigi^ hence /3g. Differences between the 
possible radius/energy relations in (|4.4I) eventually arise in the IR. However, we have verified that all 
relations in (|4.4I) lead to qualitatively similar results. For the calculations that follow, we will adopt the 
prescription (|4.2p . In this way, when N f — 0, we recover the /? function of the MN model (specifically, 
the /3 function obtained in |22)). 

It is convenient to rescale away the parameter N f in the master equation ()2.1ip by the change 
P = NfP and Q — NfQ. This leads to the following scaling for the /3 function, 

1^9 = ^^/?s(a;,3) , 9 = -y=9 , (4.7) 

where 

(4.8) 

In what follows we will thus compute (3g. Note that, in terms of the 't Hooft coupling A = g'^Nc, one has 
5^ — xX, X — Nf/Nc, and 

/3a = /(x,A). (4.9) 

12 



— ■ 05(x, q) , q = 


^'- 


^ = tanh(r) P{r) . 
9^ 





This can be compared with the NSVZ /? function (|l.ip , which in terms of A reads 

/3a^- ^ ^.^ ^ (3-x(l-7o)). (4.10) 

87r^(l- s^^J 



This agrees with the structure of the holographic /? function (|4.9|) . i.e. in the large iVc limit at fixed 
Nf/Nc it only depends on A and x = Nf/Nc. 

At this point it is useful to recall some basic facts of the NSVZ /3 function. It was suggested by Seiberg 
[27] that a conformal window for SQCD should exist for ^Nc < Nf < 3Nc, where a family of massless 
SQCD theories with Nf massless flavors develop an IR fixed point at finite coupling. All flavored gauge 
theories in the conformal window would be deconfined and chiral symmetry restored. The lower end-point 
should be considered a lower-bound on the actual value. A non-trivial IR fixed point can be found if 
X = j^ ~ 3 |24) . Indeed, using the explicit form of the one-loop anomalous dimension the vanishing of 
the /3 function requires 

--l = -lo = ^g'N, + 0{g'N!). (4.11) 

X OTT^ 

It is clear that this fixed point moves towards the strongly coupled region as x decreases from 3 to lower 
values. This assumes a small value of the anomalous dimension. As we will see below, the present 
holographic system, like the one of [5l|6] seems to involve large values of the anomalous dimension 70, in 
fact 70 = -1/2 in the UV. 

The calculations that follow use our specific choice for the embedding function S{r) = tanh (2r). 
However, the structure of the fixed points seems to be a generic property of the solutions of the master 
equation (j2.1ip with linear dilaton asymptotic and any embedding function S with S{r) — > 1 at infinity. 
This asymptotic includes previously known solutions with massless flavors. 

The linear dilaton asymptotics of these types of backgrounds preclude the emergence of an anti de 
Sitter geometry at infinity, which should be a more appropriate description near the UV fixed points. 
Despite this fact and despite the above mentioned ambiguities in the definition of the holographic beta 
function, we will find some remarkable coincidences with the expected behavior in flavored SQCD. 

4.1 Nf < 2N^ 

The /3 function for the gauge theory with massive fundamental flavors is obtained by taking the solution 
P{r) found in the previous section (see fig. [2^) and applying the formula (|4.6|) . The result is shown in 

fig. Si. 

The /3 function has a UV fixed point a.t g — 0, where it has the following behavior 

l39^-^i2N,-Nf)9'. (4.12) 

3z7r 

Remarkably, this exactly agrees with the NSVZ /3 function (jl.ip near g = 0, if 70 = —1/2 in the UV 
- where mass terms can be neglected. A similar conclusion was reached in the case of the backgrounds 
with 5 = 1 [SI [5] . This is not surprising, since in the UV our S differs from S' = 1 by exponentially 
suppressed terms, which do not affect the leading behavior in (|4.12p . It would be interesting to have an 
independent derivation of the anomalous dimension 70 by holographic methods, but presently it is not 
clear to us what the correct prescription would be|j 



''In [5] an attempt was made to compute 70 by proposing that the quartic coupling of the gauge theory should be 
identified with some quotient of the volumes of different cycles of the manifold. 
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Figure 4: ^JNfPg as a function of ^Njg, corresponding to the supergravity solutions in fig.[2^,b,c,d. (a) 

X = 7/4. (b) a; = 2. (c) x = 7/3. (d) a; = oo. 



The /? function of fig. |4^ is zero at g = 0, negative and monotonically decreasing for g > 0, thus 
implying asymptotic freedom and ordinary confinement in the IR, where g — > cx). In particular, we find 
no additional IR or UV fixed points at finite coupling. 

4.2 Nf = 2N^ 

Using the solution P{r) found in the previous section (see fig. [2}d) we determine the j3 function, shown 
in fig. lib. We can see that a non-trivial UV fixed point g = g* appears. Although we do not have the 
gravity solution that describes the missing branch g < g*, some interesting features can be inferred by 
comparing with the NSVZ /? function (|1.1[) in this UV region where mass terms can be neglected. For 
A^^ = 2Nc, the NSVZ /3 function becomes 



f3a 



^iV. 



167r2(l 



£N^ 



(1 + 270) 



(4.13) 



Again, it is consistent with our results if 70 — > —1/2 in the UV and g fiows to 5,. Moreover, since the 
perturbative NSVZ /3 function is negative near the UV fixed point at <; = 0, by continuity there must 
be at least another point g'^, with gl < g*, where the /3 function vanishes. In the simplest assumption 
that there is only one such point, this would be an IR fixed point. The resulting picture is in fact similar 
to the one proposed by Seiberg (for a discussion on the effect of mass terms see e.g. [28 ). Obviously, a 
description using massive flavors like the present one cannot describe the emergence of a conformal fixed 
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point in the infrared. However, given the presence of the UV fixed point a.t g — g,, the IR fixed point 
seems to be the simplest possibility that permits a negative beta function near g — 0. The combined 
presence of a pair of IR and UV fixed points is also a prerequisite for the existence of a mechanism in 
which the disappearance of the conformal window is due to the annihilation of a pair of fixed points [251 . 
Notice that if for Nf = 2Nc the IR fixed point appears at 't Hooft coupling A = 0(1) (as suggested 
by a naive extrapolation of (I4.1ip ). it would be very difficult to see it by means of perturbative and 
holographic techniques. 

4.3 Nf > 2Nc 

Using now the solution P{r) of fig. ^ we determine the f3 function for the case x — 7/3. This is shown 
in fig. |3t. The /? function has, like in the Nf < 2Nc case, a UV fixed point a.t g — 0, where it has the 
behavior 

^^--^^^f-^^'^^s'- (4.14) 

This exactly agrees with the NSVZ /3 function (jl.ll) of the Seiberg dual gauge theory with Nc = Nf — Nc 
near g = 0, if again we set 70 — —1/2 in the UV. This strongly suggests that in the UV region the 
background obtained with our boundary conditions describes, when Nf > 2Nc, the Seiberg dual system. 

It must be stressed that in the present case Seiberg duality is only an approximate relation that 
depends on the scale of energy (see ^Q\). The idea is that at a given scale ^ one can integrate out 
massive flavors which have mass greater than fi and remain with a reduced number of light flavors. In 
the present framework, this reduced number of flavors at an energy scale r is effectively described by 
Nf{r) = Nf S{r). With our choice of S{r), massive flavors are accumulated near r w 0.5 (see fig. [Ij. In 
the infrared region, where r ~ 0, one has 5' ^ r** so Nf{r) — > 0, as expected since in this region the energy 
scale is much smaller than the characteristic mass of the flavors. On the other hand, in the UV region, 
5 — >■ 1 and Nf (r) — >■ Nf , which is consistent with the fact that at this scale of energies all flavors look 
massless. As observed in [10_, the master equation (12.111) remains invariant under Nc — > Nf{r) — Nc and 
Nf{r) — > Nf{r). This transformation changes Q{r) — >■ —Q{r). This is the only sense in which Seiberg 
duality can be applied to the present system (in particular, Nc — )■ Nf — Nc and Nf — >■ Nf is not a 
symmetry of the master equation) and it is consistent with our proposal that the solution P{r) of flg. 
[5J; describes the Seiberg dual system at an energy scale much larger than the characteristic mass of the 
flavors, where Nf{r) — )■ Nf. 

Having obtained a gravity solution for the "Seiberg dual" system, the question is how to identify a 
background dual to the original gauge theory. When Nf > 2Nc, we expect that the gauge theory will 
develop a Landau pole. This means that the theory cannot be extended beyond a certain UV scale. On 
the gravity side, it means that the geometry should terminate at a maximum value of r, where it probably 
has a singularity. Indeed, there is a one-parameter family of solutions with parameter hi that at r = 
approach the deformed MN solution but at some flnite r meet a singularity where P — \Q\. These are 
the solutions which have an hi whose value is anything lower than the critical hi of flg. [3] In this case 
we lack a clear criterium to pick a unique solution in this family that is dual to the original gauge theory. 
It should also be noted that the application of holography is difficult to justify for singular backgrounds 
that do not get to infinity. 

The P function in fig. HJ; exhibits a local maximum precisely near the g^ where a fixed point appears 
in the Nf = 2Nc case. Indeed, as Nf approaches 2Nc, the local maximum approaches the line /3g = 
and occurs at large values of r. In the strict Nf — 2Nc limit, the branch g < g* disappears from the 
figure, because the solution gets to r = 00 already at g = g, . 
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Finally, fig. |3Ji shows the gauge coupling (3 function computed in the infinite flavor limit, that is, 
for the solution shown in fig. [2ji. It shares similar features as the case x = 7/3, except that the local 
maximum has disappeared. The disappearance of the local maximum can be understood as follows: for 
x = oo, one has x = Nj/Nc — > 1, where Nc = Nf — Nc- Thus one is computing the (3 function of the 
"Seiberg dual" system with x = 1. For cc = 1, the Pg indeed looks very similar to fig. HJi. 

5 Conclusions 

We have investigated the new gravity backgrounds found in |10| dual to A/" = 1 supersymmetric gauge 
theories with massive fundamental flavors. These backgrounds are characterized by a function S{r) which 
encodes the flavor brane distribution. In the specific backgrounds studied in this paper we have chosen 
a continuous S{r) = tanh (2r), with support in the whole space Q < r < oo, which leads to a simple 
analytic expression for the function Q{r), and thus permits a more straightforward integration of the 
master equation (|2.1ip that determines P(r), hence the complete geometry. The solutions -parametrized 
by X = Nf/Nc- were uniquely determined by imposing boundary conditions that ensure regularity at 
r = and acceptable asymptotic behavior at infinity. In this way, solutions are free from the IR singularity 
that affects the massless flavor 5=1 case of [5l [6] . 

We have then investigated properties of the gauge coupling beta function and the possible emergence 
of fixed points. As explained in [5l[6], the main feature that seems to determine the properties of the dual 
gauge theory is the presence of quartic operators in the superpotential that arise upon integration of the 
Kaluza-Klein modes of the original string theory. These operators are of the form W ~ h{Q^Qu){Q"'Qr 
with gauge indices contracted inside the parentheses and lead to a sextic potential in the scalar fields 
They become marginal when the anomalous dimension 70 is —1/2. For this value of the anomalous 
dimension the NSVZ 13 function (|l.ip becomes 



I 



3g3 i2N,-Nf) 

We have seen that this expression agrees with the holographic /3 function in the UV region for Nf < 2Nc, 
and also for Nj > 2Nc if we replace Nc — ^ Nf — N^. We argued that for Nj > 2Nc our backgrounds 
should therefore describe the Seiberg dual system, in the generalized sense discussed in sect. 4.3. After 
this replacement, the /? function stays negative for all Nf > 2Nc. As discussed, in the Nf > 2Nc case, 
finding a gravity description of the original system before Seiberg duality is difficult because, as the 
expression (j5.ip indicates, asymptotic freedom is lost and thus a Landau pole is expected, presumably 
meaning that the gravity solution must encounter a singularity at some radius r. 

For Nf — 2Nc we found that the theory has a UV fixed point, which hints at the presence of an IR 
fixed point at some lower coupling, if one is to match continuously with standard perturbative results. 
In this context, we stress that our theory only converges to the massless case in the UV, and it can only 
asymptotically recover the presence of conformal fixed points. In particular, near the would-be IR fixed 
point the theory can at most have an approximate conformal symmetry at energies greater than the flavor 
mass scale. For Nf < 2Nc we have not found any evidence of an IR fixed point, perhaps suggesting that 
in the presence of quartic operators the "conformal window" opens and closes at Nf = 2Nc. 

The effect of higher dimensional operators such as h{Q^Qu){Q^Qr) -which in the component La- 
grangian leads to terms (scalars)^ and (scalar)^ (fermion)2-has a counterpart in non-supersymmetric 



A general discussion of quartic operators can be found in |28| 
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QCD. It produces effects that are similar to well known non-perturbative effects related to chiral dynam- 
ics in the low energy effective field theory. For example, a quartic fermion operator has a key role in 
the emergence of chiral symmetry breaking and must have a role in the disappearance of conformality: 
Schwinger-Dyson gap equation for the fermion propagator implies a direct relation between the onset of 
chiral symmetry breaking, thus the presence of a non vanishing chiral condensate, and the point where 
the four-fermion operator becomes relevant in the RG flow [111 123] . This happens by lowering the number 
of flavors in QCD-like theories, starting from the point where asymptotic freedom sets in. At sufficiently 
low Nf chiral symmetry will always be broken and conformality is lost. 

An interesting open problem is finding regular backgrounds that can describe the massless flavor limit 
in a controllable manner. The current approach uses an S(r) function that determined the flavor mass 
scale r ~ Mf. One can attempt to study the limit Mf — >■ within this context. Although this approach 
seems to be affected by singularity problems similar to those of the massless S = 1 case, it is possible 
that some universal properties can be learned by studying this limit in detail. 
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